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Abstract 



Kaon production, in particular K'^ production in heavy ion collisions at intermediate energies 
I is discussed. Main emphasis is put on the question if subthreshold production can serve as a 

■ suitable tool to test the high density phase of such reactions and to deliver information on the high 

VO ■ density behavior of the nuclear equation of state. It is shown that the excitation function in 

heavy (Au + Au) over light (C + C) systems provides a robust observable which, by comparison 
to data, strongly favors a soft equation of state. A second question of interest is the existence of 
CO ■ an in-medium kaon potential as predicted by effective chiral Lagrangiens. Here it is argued that 

transport calculations support this scenario with, in the meantime, a significant level of consistency. 



g 1 Introduction 



The original motivation to study the kaon production in heavy ion reactions at intermediate energies, 
namely to extract information on the nuclear equation of state (EOS) at high densities is a matter of 
5^ ! current debate. Already in the first theoretical investigations by transport models it was noticed that 
the yield reacts sensitive on the nuclear equation of state [H |21 13 El- The yields were found to 
be about a factor 2-3 larger when a soft EOS was applied compared to a hard EOS. At that time 
the available data [S] already favored a soft equation of state. However, calculations as well as the 
experimental data were still burdened with large uncertainties. 

In jni we studied the question if in the meantime decisive information on the nuclear EOS can be 
extracted from subthreshold kaon production in heavy ion collisions. There are several reasons why it 
appears worthwhile to do this: Firstly, there has been significant progress in the recent years towards a 
more precise determination of the elementary kaon production cross sections [71 |H] , based also on new 
data points from the COSY-11 for the reactions pp — > pAK~^ very close to threshold pi. Secondly, the 
KaoS Collaboration has performed systematic measurements of the production far below threshold 
in heavy {Au + Au) and light (C + C) systems Looking at the ratios built from heavy and light 
systems possible uncertainties which might still exist in the theoretical calculations should cancel out 
to a large extent which allows to draw reliable conclusions. Furthermore, far below threshold the kaon 
production is a highly collective process and a particular sensitivity to the compression of the participant 
matter is expected. 
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2 The Model 



The present investigations are based on the Quantum Molecular Dynamics (QMD) transport model 
For the nuclear EOS we adopt soft and hard Skyrme forces corresponding to a compression modulus of 
K=200 MeV and 380 MeV, respectively, and with a momentum dependence adjusted to the empirical 
optical nucleon- nucleus potential The saturation point of nuclear matter is thereby fixed at Eb = 
-16 MeV and psat = 0.17 fm"^ [TTl- The calculations include A(1232) and A^*(1440) resonances. The 
QMD approach with Skyrme interactions is well tested, contains a controlled momentum dependence 
and provides a reliable description of the reaction dynamics in the SIS energy range, expressed e.g. by 
collective nucleon flow observables as well as particle production. In contrast to AGS energies where 
the creation of resonance matter may lead to an effective softening of the EOS, baryonic resonances 
with masses above the A^*(1440) can safely be neglected for the reaction dynamics at SIS energies. 

We further consider the influence of an in-medium kaon potential based on effective chiral models 
[121 ESI UH UH lini- The mean field consists of a repulsive vector part = 3/8/*^j^ and an 
attractive scalar part S5 = mx — = mx — \/^x ~ Sxn//^P5 + V^V^^. Here is the baryon vector 
current and ps the scalar baryon density and Sxn = 450 MeV. Following J3] in the vector field the 
pion decay constant in the medium /*^ = 0.6/^ is used. However, the enhancement of the scalar part 
using /*^ is compensated by higher order contributions in the chiral expansion J3] , and therefore here 
the bare value is used, i.e. T,knPs/ fn- Compared to other chiral approaches [121 Hlj the resulting 
kaon dispersion relation shows a relatively strong density dependence. The increase of the in-medium 
mass mx, Eq. (|2I), with this parametrisation is still consistent with the empirical knowledge of 
kaon-nucleus scattering and allows to explore in-medium effects on the production mechanism arising 
from zero temperature kaon potentials. For the kaon production via pion absorption ttB — > YK^ the 
elementary cross section of are used. For the A^A^ — > BYK^ channels we apply the cross sections 
of Ref. |Zj which give a good fit to the COSY-data close to threshold. For the case of A^A — > BYK^ 
and AA — > BYK^ reactions experimental data are rare. Thus we rely on the model calculation of 
ref. [H]. In the case that a A^* resonance is involved in the reaction we used the same cross section as 
for nucleons. In the presence of scalar and vector fields the kaon optical potential in nuclear matter has 
the same structure as the corresponding Schroedinger equivalent optical potential for nucleons 

f/opt(p, k) = + —k^v^ + . (1) 

mx 2mx 

and leads to a shift of the thresholds conditions inside the medium. To fulfill energy-momentum 
conservation the optical potential is absorbed into an newly defined effective mass 

mx(p, k) = -y/mf^ + 2mxt/opt(p, k) (2) 

which is a Lorentz scalar and sets the canonical momenta on the mass-shell = fc^ — m^. Thus, e.g., 
the threshold condition for production in baryon induced reactions reads i/i > fhB + friY + rhK with 
a/s the center-of-mass energy of the colliding baryons. For a consistent treatment of the thresholds 
the scalar and vector baryon mean fields entering into eq. ^ are determined from two versions of the 
non-linear Walecka model with K=200/380 MeV, respectively |3j. The hyperon field is thereby scaled 
by 2/3 which yields also a good description of the A flow ^HI- Since the parameterizations chosen for 
the non-linear Walecka model yield the same EOS as the Skyrme ones, the overall energy is conserved. 
The kaon production is treated perturbatively and does generally not affect the reaction dynamics |19j . 

3 Probing the nuclear EOS by subthreshold production 

The excitation function for Au + Au and C + C reactions starting from 0.8 A-GeV which is far 
below threshold (-Ethr = 1-58 GeV) has been measured by the KaoS Collaboration [1011201 • In [01 we 
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calculated this excitation function for a soft/hard EOS including the in-medium kaon potential. For 
both systems the agreement with the KaoS data [TU] is very good when a soft EOS is used. In the 
large system there was a visible EOS effect which is absent in the light system. The inclusion of the 
repulsive in-medium potential is thereby essential to reproduce the data PUI- Already in the light 
system the yield is reduced by about 50%. To extract more clear information on the nuclear EOS, 
in Fig. 1 we considered the ratio R of the kaon multiplicities obtained in Au -t- Au over C + C reactions, 
normalized to the corresponding mass numbers. The kaon potential is included since without the in- 
medium potential one is not able to reproduce the experimental yields jSl EI] • The calculations are 
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Figure 1: Excitation function of the ra- 
tio R of K'^ multiplicities obtained in 
inclusive Au + Au over C + C reac- 
tions. The calculations are performed 
with in-medium kaon potential and us- 
ing a hard/soft nuclear EOS and com- 
pared to the data from the KaoS Col- 
laboration pn] . 



performed under minimal bias conditions with 6max = 11 fm for Au + Au and femax = 5 fm for C + C and 
normalized to the experimental reaction cross sections fOll^ni- Both calculations show an increase of 
R with decreasing incident energy down to 1.0 A- GeV. However, this increase is much less pronounced 
when the stiff EOS is employed. In the latter case R even decreases at 0.8 A-GeV whereas the soft 
EOS leads to an unrelieved increase of R. At 1.5 A- GeV which is already very close to threshold the 
differences between the two models become small. The strong increase of R can be directly related to 
higher compressible nuclear matter. The comparison to the experimental data from KaoS ^U] where 
the increase of R is even more pronounced strongly favors a soft equation of state. 

To obtain a quantitative picture of the explored density effects in Fig. 2 the baryon densities are 
shown at which the kaons are created. The energy is chosen most below threshold, i.e. at 0.8 A-GeV 
and only central collisions are considered where the effects are maximal. dM^+ldp is defined as 

'^^'+ dP- 

dMK./dp = E , , \, (3) 

where is the baryon density at which the kaon i was created and Pj is the corresponding production 
probability. For the comparison of the two systems the curves are normalized to the corresponding mass 
numbers. Fig. 2 illustrates several features: Only in the case of a soft EOS the mean densities at which 
kaons are created differ significantly for the two different reaction systems, i.e. < p/psat >=1.46/1.40 
for C + C and 1.47/1.57 for Au + Au using the hard/soft EOS. Generally, in C + C reactions densities 
above 2psat are rarely reached whereas in Au + Au the kaons are created at densities up to three times 
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Figure 2: Kaon multiplicities (normal- 
ized to the mass numbers of the collid- 
ing nuclei) as a function of the baryon 
density at the space-time coordinates 
where the mesons have been cre- 
ated. Central (b=0 fm) Au + Au and 
C + C reactions at 0.8 A-GeV are 
considered. The calculations are per- 
formed with in-medium kaon potential 
and using a hard/soft nuclear EOS. 



saturation density. Furthermore, for C-l-C the density distributions are weakly dependent on the nuclear 
EOS. The situation changes completely in Au + Au. Here the densities profile shows a pronounced EOS 
dependence Moreover, the excess of kaons obtained with the soft EOS originates almost exclusively 
from high density matter which demonstrates that compression effects are probed. 
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Figure 3: Kaon multiplicities func- 
tion of Apart at E=1.0 AGeV. Multiphc- 
ities as well as Apart are normalized to 
the mass numbers of the colliding nu- 
clei. 



The density effect is also clearly reflected in the Apart dependence of the kaon multiplicities. Fig. El 
compares the Apart dependence in Au -f Au and C + C reactions, now for 1.0 AGeV laboratory energy. 
In both cases the kaon multiplicities as well as Apart which has been derived within the geometrical 
model, are normalized to the corresponding mass numbers. As can be seen from there, in C + C the 
Apart dependence of the kaon production is completely insensitive to the nuclear equation of state. The 
large system, in contrast, shows a distinct EOS dependence. In Au + Au the enhanced kaon production 
which is due to higher compression when a soft EOS is used, becomes more and more pronounced with 
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increasing centrality. However, one has to keep in mind that the large difference between the soft and 
hard EOSs in most central reactions is washed out to some extent in minimal bias reactions. There 
the bulk of kaons originates from semi-central reactions 6 ~ 5 fm, corresponding to Apart/^max ~ 0.7. 
In this context it will be very helpful to study the EOS dependence also in C + Au reactions and to 
compare to forthcoming data from KaoS 

3.1 How firm are the conclusions? 

Of course now the question arises, how firm conclusions on the nuclear EOS are which can be drawn 
from kaon production in heavy ion reactions. Possible concerns might be based on the facts that 
subthreshold kaons are an extremely rare probe and not all elementary production cross sections have 
been measured. 

In heavy ion collisions the K'^ production runs over two major channels, namely baryon-baryon 
induced reactions BB i— BYK^ and pion-baryon induced reactions ttB YK^ which are both about 
equally important (23- In both cases the initial baryons can either be nucleons or nucleon resonances 
(mainly A(1232)), the hyperons are A or S hyperons. Processes with nucleon resonances in the final 
state are energetically suppressed. Concerning the knowledge of the elementary reaction cross section 
the situation is presently as follows: the NN and ttN cross sections are quite well under control since 
these channels have been measured in pp reactions jU] and in vr^p reactions. The reactions which involve 
nucleon resonances, in particular with A's in the initial states {i = NA, vrA, AA) are less secure due 
to the lack of corresponding experimental data. Thus one has to rely on model assumptions. The cross 
sections which have been used in the present transport calculations are based on the effective Lagrangian 
model of Refs. [HI HZl |21j . The isospin dependence of the cross sections is thereby determined in the 
standard way by isotopic relations assuming iso-spin independent matrix elements. 

Hence there exists still some uncertainty in the transport calculations due to the incomplete knowl- 
edge of the elementary reaction cross sections. There are, however, two good arguments why conclusions 
on the EOS dependence of the kaon production should be rather robust against such possible uncer- 
tainties: 

• Changes of the production cross sections shift absolute yields but considering the ratio of different 
reaction systems such errors drop out in leading order. 

• Conclusions are based on the slope of this ratio as a function of energy. It is rather unlikely 
that an incomplete knowledge of the cross sections, e.g. concerning their isospin dependence, can 
create the observed energy dependence. The systematics of spurious contributions should be flat 
as a function of energy. Otherwise one have to assume extremely unconventional threshold effects. 

To be more quantitative we consider in the next figure the excitation function of R for the various 
production channels. There the ratios Ri are built separately for the production channels with initial 
states i = NN, nN, NA, vrA, AA. As can be seen from Fig. |3J the shape of R is not strongly influenced 
by the A^A , vrA channels which are the most insecure ones. The excitation function for the A^A 
contribution varies only little as a function of energy and is similar using the different EOSs. The 
contribution of the ttA channel is decreasing for both, a hard and a soft EOS. The shape of R is 
to most extent determined by the A^A^ and vrA^ contributions. In our calculations the latter channel 
is responsible for the decrease of R very far below threshold when the hard EOS is applied. These 
findings are generally confirmed by independent transport calculations of the Nantes group using the 
IQMD transport model [21] shwon in Fig. 5. These calculations include an in-medium kaon potential 
derived in relativistic mean field theory (RMF) |25] which is somewhat less repulsive than that one 
used in our calculations. For the soft EOS the IQMD calculations coincide almost with the present 
results 0. For the hard EOS there exist still deviations concerning the slope of R going far below 
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Figure 4: Dependence of the excitation function of R on the various production channels. Central 
(b=0 fm) Au + Au and C + C reactions are considered. The calculations are performed with in-medium 
kaon potential. 

threshold. This could be due to the different in-medium potentials and is an open question which has 
to be resolved by future investigations. However, the two sets of transport calculations show a good 
overall agreement and both rule out the hard EOS from the comparison with data. The shaded area 
in Fig. 5 can be taken as the existing range of uncertainty in the theoretical model description of the 
considered observable. 

Moreover, the IQMD calculations were also repeated with an alternative set of A^A; A A i— s> NY 
cross sections taken from [21] which are about a factor of two larger than those from Tsushima et al. 
[S]. The ratio R is almost completely independent on this change. Another prove for the robustness of 
this observable is the fact influence of the repulsive potential which decreases the total kaon yield 
by about a factor drops almost completely out when ratio of the two different mass systems is built 

urn- 

3.2 Are the conclusions consistent with information from other sources? 

Finally the question arises up to which degree a consistent picture has emerged after more than ten 
years of intensive experimental and theoretical efforts to understand the kaon production at intermediate 
energies. In the following I will argue that concerning we have in the meantime a rather consistent 
picture while for K~ the situation is not yet so clear. 

The reason is that subthreshold production is easier to handle, both from the experimental and 
the theoretical side. Due to the lower threshold there are much more data with much higher precision 
available than for K~ . Also theoretically the medium dependence of the meson properties are better 
under control. The mean field approximation seems to work, i.e. mass shifts and can be taken from 
the leading order chiral Lagrangian [T^ IT^ and in the medium there exits still a well established quasi- 
particle pole [23 El- This allows a treatment within standard transport based on the quasi-particle 
approximation. The K~ — N system, on the other hand, lies in the vicinity of the A1405 resonance which 
implies a strong coupling to this state. Hence a simple mean field picture will not work but a coupled 
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Figure 5: Excitation function of the 
ratio R of multiplicities obtained 
in inclusive Au + Au over C + C re- 
actions. Our results are compared to 
independent IQMD calculations pij . 
The shaded area indicates thereby the 
range of uncertainty in the theoreti- 
cal models. In addition IQMD results 
based on an alternative set of elemen- 
tary production cross sections are 
shown. 
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channel treatment of the K~ in the medium is necessary |27| |2Hl I2S]- As a consequence, the K~ has 
complicated spectral properties in the medium and a simple quasi-particle picture is nor more suitable 
)28L ■ This requires a more sophisticated treatment within transport simulations which accounts at 
least approximately for the off-shell contributions of the spectral functions [30j. In the latter case the 
interpretation of existing data with the help of transport simulation has not yet reached a level which 
allows to say that the properties are settled. 

Coming back to the the situation is much more satisfying. Most transport calculations agree 
on the necessity of a repulsive potential in order to understand total yields as well as the collec- 
tive motion of mesons, i.e. in-plane and out of-plane flow 311 ISH ISSl OH ISHl ES|- This picture 
was recently complemented by measurements of the production in proton-nucleus reactions |37j . 
Although such reactions test only subnormal nuclear densities they are much easier to handle than the 
complicated dynamical evolution of heavy ion reactions. Corresponding data pA from ANKE revealed 
strong evidence for a repulsive potential which is of the order of magnitude as predicted by effective 
chiral lagrangiens. 

3.2.1 Many-body calculations 

Concerning the nuclear equation of state one has to confront the information from subthreshold 
production with the knowledge obtained from other sources: At intermediate energies heavy ion re- 
actions test the density range between two and three times nuclear density. The information from 
kaon production implies that in this density range the EOS should show a soft behavior. One has of 
course to be aware that the adopted Skyrme forces are simplified interactions which are easy to handle 
but must not be very realistic. A microscopic approach to nuclear matter which is based on realistic 
nucleon-nucleon interactions has to solve the correlated quantum-mechanical many-body problem. Such 
an approach is ,e.g., the Brueckner-Hartree-Fock approach which accounts for the two-body correlation 
in the medium solving the Bethe-Goldstone equation, i.e. the Lippmann-Schwinger equation in the 
medium. The relativistic version, i.e. the Dirac-Brueckner-Hartree-Fock (DBHF) approach delivers 
quite reasonable values for the nuclear saturation mechanism jSH] and can thus be considered as a reli- 
able method to extrapolate the many-body approach to higher nuclear densities. In Fig. 6 we compare 
DBHF results obtained with two slightly different NN-interactions, i.e. Bonn A and Bonn B |39j to the 
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Figure 6: The Nuclear mat- 
ter EOS from microscopic Dirac- 
Brueckner-Hartree-Fock calcula- 
tions from 110] are compared 
to that resulting from soft /hard 
Skyrme forces. 



-20 







1 



2 



3 



4 







simple Skyrme type equations of state. The results of advanced DBHF calculations are taken from |4()j . 
In this context it may worthwhile to mention that very similar results to [lOj were recently obtained 
by Weise and coworkers in a treatment of the nuclear many-body problem based on chiral perturbation 
theory pi]. Generally such many-body calculations predict a relatively soft behavior of the EOS in the 
relevant density range, i.e. for densities below about three times saturation density. In the microscopic 
approach the high density behavior is, however, only loosely connected to the curvature at saturation 
density. For Bonn A and B the compressibilities are quite different while the EOSs at high densities 
are very close. Below 3po both are not too far from the soft Skyrme EOS. 

3.2.2 Nucleon Flow 

Another observable which helps to constrain the nuclear mean field and the underlying EOS at supra- 
normal densities is the collective nucleon flow The transverse flow vi has been found to be sensitive 
to the EOS and, in particular in peripheral reactions, to the momentum dependence of the mean fleld 
I44j . The elliptic flow V2, in addition, is very sensitive to the maximal compression reached in the 
early phase of a heavy ion reaction. The cross over from preferential in-plane flow f 2 < to preferential 
out-off-plane flow f 2 > around 4-6 AGeV has also led to speculations about a phase transition in this 
energy region which goes along with a softening of the EOS 

The present situation can be summarized as follows: At SIS energies existing flow data are consistent 
with the usage of mean flelds which are close to those obtained from microscopic DBHF calculations, 
both concerning their density and momentum dependence E^ . A detailed comparison to FOPI 
data 46j for vi and V2 between 0.2 and 0.8 AGeV favors thereby a relatively soft EOS [44j such as the 
DBHF result for Bonn A (K=230 MeV, shown in Figure 6). The full flow excitation function, ranging 
from low SIS up to top AGS energies, has been studied in The conclusion from this study was 
that, both, super-soft equations of state (K=167 MeV) as well as hard EOSs (K>300 MeV) are ruled 
out by data. Hence the picture is again consistent with the information obtained from kaon production. 

4 Summary 

To summarize, we flnd that at incident energies far below the free threshold K"^ production is a suitable 
tool to study the dependence on the nuclear equation of state. Using a light system as reference 
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frame there is a visible sensitivity on the EOS when ratios of heavy {Au + Au) over hght {C + C) 
systems are considered. Transport calculations indicate that the production gets hardly affected 
by compressional effects in C + C but is highly sensitive to the high density matter (1 < p/psat < 3) 
created in Au + Au reactions. Results for the K'^ excitation function in Au + Au over C+C reactions as 
measured by the KaoS Collaboration, strongly support the scenario of a soft EOS. This statement turns 
out to be rather robust against possible model uncertainties: It is almost independent on the variation 
of particular reaction channels where elementary cross sections are uncertain. It is also insensitive to 
the inclusion/neglection of a changing in- medium kaon mass as predicted by chiral models. The idea 
of a soft EOS in the considered density range is also consistent with the knowledge from microscopic 
many-body theory and from nuclear flow analysis in heavy ion reactions. 

Concerning the quest for in-medium modification of the kaon properties transport calculation have 
in the meantime reached a certain level of consistency: The explanation of the total yields and the 
K'^ flow requires the presence of a repulsive in-medium potential. This picture has been complemented 
by measurements of kaon production in p + A reactions. 
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